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Abstract 

An output feedback sliding mode controller (SMC) is proposed to minimize the effects of vibrations of 
slewing flexible manipulators. A spline trajectory is used to generate ideal position and velocity commands. 
Constrained nonlinear optimization techniques are used to both calibrate nonlinear models and determine 
optimized gains to produce a rest-to-rest, residual vibration-free maneuver. Vibration-free maneuvers are 
import ant for current and future NASA space missions. This stud y required the development of the non- 
linear dynamic system equations of motion; robust control law design; numerical implementation; system 
identification; and verification using the Sandia National Laboratories flexible robot testbed. Results are 
shown for a slewing flexible beam. 


Introduction 

For NASA space applications, lightweight robotic manipulators are necessary to reduce launch costs, 
power consumption, and storage volume of the robot. Slewing structures with long flexible members, such as 
the Shuttle Remote Manipulator System (RMS) and the Space Station RMS, can excite vibrations. These 
vibrations can severely degrade the pointing accuracy, thus limiting the speed of rotation and productive 
use of current and future robotic systems. To achieve good control performance and position precision with 
current technology requires massive stiff manipulators. Since mass is the strongest driver of launch costs, 
massive telerobotic systems are unacceptable. Inherent flexibility for manipulator systems is a consequence 
of launch mass minimization. Flexibility is difficult to model but without its inclusion in the dynamic 
model, slewing performance will remain poor and cent rol marginally stable. The focus of this research is 
the development of a robust control system that demonstrates residual -vibration suppression and robust 
tracking using only colocated joint sensors and actuators. The methodology includes development of the 
dynamic system equations of motion; sliding mode control system design; optimized model matching and 
gain calculation; and experimental verification using the Sandia National Laboratories flexible robot testbed. 

To find a practical feedback control for flexible arms, many researchers have investigated various control 
methods. A review of some of these techniques is given in Yeung and Chen . 1 Yeung and Chen also demon- 
strated successful feedback control of flexible arms using the sliding-mode technique. Nathan and Singh 2 
developed a design approach for the control of a flexible robotic arm using variable structure system theory 
and pole assignment technique for stabilization. The closed-loop system was robust to variations in payload. 
Qian and Ma 3 have introduced variable structure sliding-mode technique for tip position control. The con- 
troller performance was evaluated through simulations. Choi and Shin 4 developed a sliding mode controller 
for tip position control of a single-flexible link manipulator subjected to parameter variation. Their algo- 
rithm showed fast and favorable system responses while maintaining low sensitivity y to imposed uncertainties. 
Parker and Robinett 5 developed an output feedback sliding mode control approach for nonlinear systems in 
general, with applications to flexible manipulators. Asymptotically stable sliding surfaces are specified in the 
output space. In addition, a constraint was derived, based on L yapunov’s direct method, ensuring stability 
of the closed-loop system. The results presented in this paper is a realization of the output feedback sliding 
mode controller. 
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Dynamic Model 


The dynamic equations of motion for both the rigid body, 6 and the flexible body, q i (t) degrees-of-freedom 
(DOF) are found using Lagrange’s equations; 
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The Lagrangian is, L-T- V + Wf, where the kinetic energy, T, the potential energy, V, and the work 
from external forces, Wf are defined as 
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where i and j = 1 are the number of flexible DOF and x is the location along the beam. Applying Lagrange’s 
equations results in the following nonlinear equations of motion; 


M(x)x + N(x, x) + K(x, x) = B(x)U ( 6 ) 


y = Cx ( 7 ) 

where, x is an nx 1 vector of total DOF’s, M is a nxn configuration dependant mass matrix; N is a nx 1 
vector of Coriolis and centripetal acceleration terms; K is a rum configuration dependant stiffness matrix 
including centrifugal stiffening terms; B is a nxm matrix of control weighting coefficients; U is a mz 1 vector 
of torque inputs; y is an rarl vector of measurable outputs, C is an rxn matrix relating state variables to 
measurable outputs. 

The dynamic equations of motion were developed using the method of quadratic modes. 6 Figure 1 shows 
a schematic of the slewing flexible beam defining the mathematical geometry. 



Figure 1: Slewing Flexible Beam Schematic 
An expression for the deformation of a point along the beam is 

u(x,t) =u x (x,t)bi +y(x,t)h (8) 
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where fe f are unit vectors associated with a moving coordinate system attached to the hub. The n»- unit 
vectors are associated with an inertial coordinate system. Define the following relationships for axial and 
transverse deflections as 


u x (x,t) - 
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y(x,t)=^2^(x)q l (t) . 
1 = 1 
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Equation (10) is assumed to be separable into <j>'{x), the mode shape basis functions and q'(t), the corre- 
sponding time-dependent generalized coordinates. The following equation gives the velocity of each point 
along the rotating beam’s length: 


N d ( 1 

v(x, t ) — {[ r + X ] b l +«(*! t)j ■ 


( 11 ) 


Performing the mathematical expansions and substituting the expressions for the kinetic energy, strain 
energy, and external work into Lagrange’s equations (1,2), we arrive at the following equations for beam 
deflection and rotation, respectively. 


El [ (x)&" (x) dx 

Ja 



<j> i {x)<td(x)dx + 2p f [r+ x]g ij {x) dx 6 2 q* 
jOi 


} 


+ 



q' + 



r L . 

= T I ^(x^dx 
Jo 


( 12 ) 


+ prL 2 + pr 2 L + M t (r + L)^0 +j/> J (r + x)<f>'(x) dx q 1 = r (13) 

where p =p +M t 6(x- L) and p is mass per unit length. 

Cantilever mode shapes given by Blevins 7 , were used for this analysis and the quadratic modes for a beam 
were defined as; 

= -\f v'izw'tm- (i4) 

The final equations were arranged into the form of equations (6) and (7). 


Sliding Mode Control 


Sliding mode control provides an alternative to robot control with unknown parameters. The main advan- 
tage of SMC is its robustness to input disturbances once the sliding surface is reached. SMC uses a strategy 
whereby the active control law at any given ti me is chosen from a predefi ne set of control I aws based on 
the current state of the system. SMC takes advantage of control law switching to move a system from an 
initial state to a prescribed surface in the state space. Once on that surface, a second control law is used 
to keep the state from leaving the surface while moving toward the desired final state. Using Lyapunov’s 
direct method, SMC has been shown to be stable 5 . Furthermore, it is robust to model-parameter uncertainty 
and disturbances if bounds are known a priori. Sliding surfaces are identified for each sensor leading to the 
output feedback control law. A thorough development of SMC including several practical examples can be 
found in Utkin . 8 

The sliding surface may be chosen as 


s = W(y — y r ) + (y — y r ) = 0 (15) 

where y r is the desired sensor output time history and W is a positive definite matrix with real valued 
elements. The equivalent control is found by enforcing a condition of stationarity on the sliding surface, 
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s = W(y - y r ) + (y -y r ) = 0. (16) 

substituting for y from the equations of motion equations (6) and (7), intos and solving for CM(x)- ^(xjU 
yields; 


CM(x) -1 B(x)U = CM(x) -1 N(x, x) + CM(x) -1 K(x, x) + y r - W(y- y r ) (17) 

All terms involving x will be approximated with x, such that 

x= C*y (18) 

where C* takes on a form of a psuedo-inverse. I n addition to substituting for the esti mate of x the term 
A tanh'^ifis) is added to drive the output to the stable sliding surface ofs and the hyperbolic arctan is used 
to eliminate chatter through a boundary layer whose slope can be adjusted with /3. This results in the final 
output feedback sliding mode controller; 


u =[CM(x)' 1 B(x)]' 1 [CM(x)‘ 1 N(x,i) + CM(x) 1 K(x,x) + y r - W(y - y r ) - Atanh ^/^.(lg) 

Inversion is ensured by setting m =r. Stability has been established by using Lyapunov’s direct method . 5 
During the actual implementation of the SMC algorithm the following term was set to zero; 

CM(x) -1 N(x,i) +CM(x) -1 K(x,i) =0. (20) 

Optimization 


A constrained optimization problem was formulated for the slewing flexible beam involving the physical 
parameters of the previously derived model and an experimental response to the trajectory input. Solv- 
ing the trajectory optimization problem involved the use of a recursive quadratic programming alogrithm 
implemented in the MATLAB optimization toolbox. 9 A cost function of the form 


J = 



• Q\dt + W 2 



• e2 dt 


( 21 ) 


subject to a number of inequality constraints G(x) < 0 was used for both the model matching and the 
optimized gain analysis. 

The model matching optimization errors were set up as ej =(0 model - $ exp ) and e 2 ={eroot modcl -troot eicp ), 
where t roo t is the strain at the root of the beam. The first step involved setting the parameters belonging 
to the mass and stiffness properties of the hub and flexible link, The weights wereset to WI- 1.0 and 
W 2 = 0.0. The optimizer was allowed to formulate error predominantly during the rise time to = 0.0 to tj = 
0.35 seconds. After sufficient iterations the parameters would converge to nominal values. The second step 
was then to concentrate on the friction coefficients as parameters. The optimizer was set-up to work over 
the settling time portion of the trajectory from to =0.35 tof; =0.6 seconds until the parameters would 
converge to a nominal value. The third step used a representative set of parameters from both steps one 
and two. This set of parameters was allowed to only vary between ±20 percent. The optimizer worked over 
the combined time range of the trajectory from t 0 =0.0 to tj =0,6 seconds. Upon successful convergence 
of these parameters, the rigid body portion of the single flexible link system is identified. The final step 
included setting W 2 > 1.0 and investigating parameters directly associated with the strain location and 
beam coefficients. This resulted in closer agreement to the experimental setup but was considered only a 
second-order effect. 


For the controller optimization the errors are specified as ei =(0 C - 6) and e 2 =(0.0 - e root ) where both 
variables #and e r oot are from the simulation model. The optimizer was set up for the hub angle from to = 0.0 
totj =0.6 seconds and for the root strain from to =0.4 to tj =0.6 seconds. The cost associated with the 
root strain is for after the maneuver is completed to minimize residual vibration. For all runs WI = 1.0 and 
W 2 = 10. To start out, large steps were taken to identify possible mini mums. Starting with these mini mums 
the step size was reduced until convergence. These gains were then implemented on the hardware to obtain 
experiment al responses, 
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Experimental and Numerical Results 


TheSandia National Laboratories flexible robot testbed consists of modular flexible link/motor/hub 
mounting assemblies; elect ric DC motors and amplifiers; incremental encoders; bending strain gauges; and 
adSPACE 10 real-time control computer and data acquisition system. The slewing flexible beam parameters 
are given in Table 1. 


Parameter Symbol 

Value 

Unit 

Length 

L 

48.42 

cm 

Width 

w 

7.62 

cm 

Thickness 

t 

0.1574 

cm. 

Hub Radius 

r 

8.89 

cm 

Mass Density 

Pm 

2700 

kg/m 3 

Ti p M ass 

M t 

0,0 

kg 

Beam Stiff 

El 

0.176 

kg . m 2 

Motor Inertia 

J m 

6. 92 3 

kg . m 2 

Viscous Damp 

Kf 

1 . 3 7 4 

kg . m 2 /s 


Table 1: Slewing Flexible Beam Physical Parameters 

A numerical simulation was developed that realized the mathematical models developed earlier. MATLAB 
was used to implement the differential equations. The dynamics of the plant were treated as continuous states, 
while the control laws were treated as discrete states. All sampling was performed at 1000 Hz. 

The reference motion trajectory of the hub is generated from a spline fit of the initial hub angle, equal to 
-90° to the final hub angle, equal to 90°, for the single flexible link case. The time for all the trajectory 
runs was specified as AT = 0.35 seconds. 

Model Calibration/Matching 

The goal of this section was to identify a model that best captured the dynamics of the actual system. 
A simple PD controller was used to slew the beam. By following the steps outlined in an earlier section, 
the following plots show the match between the model and the experimental set-up. Using empirically 
determined gains Figures 2 and 3 show the calibration plots for hub angle, hub velocity, root strain and 
mid-span strain, respectively. 


Model Matching Response Model Matchin 9 Res P° nse 



Figure 2: Hub Angle and Velocity Calibration Results 
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Model Matching Response Model Matching Response 



Figure 3 : Root and Mid-Span Strain Calibration Results 
Optimized Gains for Controllers 

The calibrated model was used to predict the performance of the experimental set-up by using the gains 
determined from the constrained nonlinear optimization design. The results for the sliding mode control, 
where the W and A gains were optimized are shown in Figures 4 and 5, for the hub angle, hub velocity, 
root strain and mid-span strain, respectively. 


Optimized Gains Sliding Mode Control Response 



Figure 4: Flub Angle and Velocity Optimized SMC Results 


Conclusions 

A sliding mode controller was successfully demonstrated to minimize the effects of vibrations of slewing 
flexible beams. Optimization techniques were successful employed to determine meaningful nonlinear time 
domain models and optimized gain determination. In turn the optimized gains were used to predict flexible 
beam performance during large angle slews. These optimized gains were experimentally verifedon theSandia 
National Laboratories flexible robot testbed. The SMC architecture showed minimum residual-vibration 
suppression and robust tracking using only colocated joint sensors and actuators. Future work will involve 
the use of 1) piezoceramic strain sensors and actuators 11 ’ 12 to enhance stability and tracking performance, 
and 2) the use of a two DOF planar flexible manipulator. 
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